Neel transition of lattice fermions in a harmonic trap: a real-space DMFT study 
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We study the magnetic ordering transition for a system of harmonically trapped ultracold fermions 
with repulsive interactions in a cubic optical lattice, within a real-space extension of dynamical 
mean-field theory (DMFT). Using a quantum Monte Carlo impurity solver, we establish that an- 
tiferromagnetic correlations are signaled, at strong coupling, by an enhanced double occupancy. 
This signature is directly accessible experimentally and should be observable well above the critical 
temperature for long-range order. Dimensional aspects appear less relevant than naively expected. 

PACS numbers: 67.85.-d, 03.75.Ss, 71.10.Fd, 75.10.-b 



Starting with the achievement of Bose Einstein con- 
densation, ultracold atomic gases have led to fascinating 
insights in quantum many-body phenomena With 
the recent experimental successes in realizing quantum 
degenerate [1] and strongly interacting Q Fermi gases 
in optical lattices, such systems are widely considered as 
highly tunable quantum simulators of condensed matter 
Q : in principle, most solid-state Hamiltonians of interest 
are accessible by selecting appropriate fermionic flavors, 
interactions, and lattice parameters @. The recent ob- 
servation of the fermionic Mott transition in binary mix- 
tures of ''"A' in three-dimensional simple-cubic lattices 
[^,13] marks important progress in this respect. 

A major current challenge is the realization and detec- 
tion of quantum magnetism, i.e. spontaneous magnetic 
order, in ultracold atoms. For bosons in optical lattices, 
correlated particle tunneling and superexchange 3 - 
basic mechanisms underlying quantum magnetism - have 
already been observed. Moreover, indications for a ferro- 
magnetic Stoner instability have recently been found in 
a trapped spin-1/2 Fermi gas @. However, the antiferro- 
magnetic (AF) phase, which is generic for Hubbard-type 
models and many classes of solids at low temperatures, 
has not yet been seen in cold fermions in optical lattices. 
So far, experimental progress in this direction is mainly 
sought by trying to cool the systems down to a central 
entropy per particle s < log (2)/ 2 required for long-range 
AF order on a cubic lattice [lOl . Ill | (which is a factor of 
2 below current average entropies T^] ) ■ 

Due to the intrinsic inhomogeneity of trapped atomic 
clouds, even the qualitative interpretation of experimen- 
tal data may rely on corresponding quantitative simu- 
lations (e.g. of cloud compressibilities across the Mott 
transition [^). An appropriate quantitative theory of AF 
ordering in fermions in an optical lattice should there- 
fore capture strong correlation effects as well as the spa- 
tial inhomogeneity at the temperatures of experimen- 
tal interest. The dynamical mean-field theory (DMFT) 
is well established powerful, nonperturbative ap- 

proach to interacting Fermi systems in three dimensions 
|13llT^ . Spatial inhomogeneities of the optical lattice can 
be captured within a real-space extension of the method 



(RDMFT) [ll,[Ti| or, approximately, by applying DMFT 
within a local density approximation (LDA). In either 
case, the accessible parameter ranges and the numerical 
accuracy depend on the method chosen as DMFT impu- 
rity solver. So far mainly the numerical renormalization 
group (NRG) has been used in this context, which is 
reliable at low temperature T, but leads to artifacts at 
elevated T Thus, NRG is only adequate for studying 
effects with weak temperature dependence, such as the 
paramagnetic Mott transition [3, [l3l : or for ground state 
investigations [T6j . In contrast, quantum Monte Carlo 
(QMC) based methods are precise and even numerically 
cheap at or above the T ranges relevant for AF ordering 
and accessible in cold-atom experiments. 

Observables routinely measured in solids, e.g., using 
x-ray or neutron diffraction, may be hard to access in 
cold-atom based systems. In particular, the detection 
of the AF order parameter, although in principle possi- 
ble via noise correlations 18 1 or Bragg spectroscopy [l9j . 
is highly nontrivial. It is, therefore, essential to iden- 
tify fingerprints of AF phases that are easily accessible 
in current experiments. Ideally, such observables should 
also be sensitive to precursor effects, monitoring the ap- 
proach to the ordered phase. 

In this Letter, we present Hirsch-Fyc QMC [l^] based 
RDMFT studies of trapped spin-1/2 lattice fermions at 
low to intermediate T, employing a massively parallel 
code [2l|. We demonstrate that the onset of AF cor- 
relations at low T is signaled, for sufficiently strong 
interactions, by a significantly enhanced double occu- 
pancy D. This signature, which had been missed so 
far even in studies of the homogeneous Hubbard model, 
appears ideal in the cold-atom context, since (i) D is 
directly accessible experimentally and (ii) its enhance- 
ment measures nearest-neighbor AF correlations which 
are expected as a precursor effect already above the Neel 
temperature. We find strong proximity effects (beyond 
LDA) at the interfaces between AF and paramagnetic 
regions, implying that RDMFT is important for quanti- 
tative studies of ordering phenomena in inhomogeneous 
systems. Our detailed predictions provide essential guid- 
ance to experimentalists. 



Model and Methods - Balanced mixtures of spin-1/2 
fermions in an optical lattice are well described by the 
Hubbard model (with trapping potential Vi — Vor|/a^), 

{ij),cr i i,<7 

Here, fii^ = cl^c^^, c^^ {^1^) s-i'C annihilation (creation) 
operators for a fermion with (pseudo) spin a G {t, ].} 
at site i (with coordinates r^), t is the hopping ampli- 
tude between nearest-neighbor sites (ij), C/ > is the 
on-site interaction, and fi is the chemical potential. In 
the following, we choose Vq = 0.05t for a nearly realistic 
confinement j22| and use t as the energy unit. 

The RDMFT approach [H, ^ is based on the follow- 
ing expression for the Green function matrix G for spin 
a and Matsubara frequency iujn = *(2n -|- 1)t:T: 

[GaiiuJn)] = Uj + [iujn + n~Vi~ I]jcr(ja;„)] % . (2) 

Here the only approximation is the DMFT assump- 
tion of a local (i.e. site-diagonal) self-energy E, which 
corresponds to the momentum independence of S in 
translation-invariant systems [l^. The standard DMFT 
impurity problems, one for each inequivalent lattice site, 
are solved using the Hirsch-Fye QMC algorithm [20| : 
Trotter errors in the raw data are eliminated by ex- 
trapolation of the imaginary-time discretization At — > 
d^lil]. Due to the relatively high T considered here, sta- 
tistical errors can be avoided as well in full enumerations 
of the Hubbard-Stratonovich fields, which improves the 
data for ordered phases at large U > 12t. The resulting 
precision is much better than achievable with NRG [2l|. 
This remains true within the "slab" approximation, in 
which the full system is reconstructed from the properties 
of a central set of planes with periodic boundary condi- 
tions in all directions [2lj. The entropy S is not directly 
accessible in QMC simulations, therefore we use the ther- 
modynamic relation S{^i,T) = T^dfi' dN{^i' ,T)/dT 
(for total particle number N) 

Results ~ We consider a fermionic cloud confined in a 
cubic lattice for i^i = U /2, i.e. for half filling at the center 
(at least within LDA) and focus, initially, on a moder- 
atly strong interaction U = 12t, equal to the bandwidth 
W. In Fig. [U the core shows a nearly perfect staggered 
magnetization (first row) at low T (left column). With 
increasing T (from left to right), both the extent of the or- 
dered region and its polarization decrease until the order 
is lost at the bulk Neel temperature Tn « OAdt. Unfor- 
tunately, this most obvious signature of AF order is not 
directly accessible experimentally, since single-site resolu- 
tion [2^ has not yet been achieved for three-dimensional 
systems. The Neel transition is also not detectable via 
the particle density profiles (third row), as they hardly 
change at this scale. In contrast, the double occupancy 
Di = {ni^riii) provides a pronounced signal: at high T, 
it is featureless in the center, with a maximum value of 
about 0.025. Only at low T, it is enhanced, by up to 50%, 
in the emerging central antiferromagnetic core. This sig- 
nal remains strong even after integrating over the line of 
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FIG. 1; (Color online) First row: AF order in central plane of 
atomic cloud (in cubic optical lattice) for intermediate cou- 
pling U = W = 12t. At low T < 0.2t (left column), a 
large AF core is strongly polarized (nearly 90%); with in- 
creasing r, both extent and magnitude of the AF order de- 
cay. Second row: the double occupancy (upper/lower half: 
central/column density) is strongly enhanced in AF regions; 
in contrast, the density (3'^'* row, central/column) shows little 
T dependence. 4**^ row: entropy densities (central/column). 

sight (z axis) as in most experimental detection schemes; 
resulting column density plots are shown in the lower 
halves of the split images in Fig. [TJ We note that total 
double occupancies have already been measured with rel- 
ative accuracies of 20% Q ; a higher precision and spatial 
resolution appear feasible. Entropy, on the other hand, 
cannot be measured in situ; thus, its suppression in the 
AF core (Fig. [1] bottom) is not visible experimentally. 

The observed enhancement of D is also present in the 
homogeneous Hubbard model at half filling (cf. Fig. [2]), 
for which the underlying mechanism is most easily ex- 
plained. For high T thermal fluctuations prevail, which 
increase D with increasing temperature towards the un- 
correlated limit Di — >■ {ni^){nii); this effect has already 
been used in thermometry (comparing, e.g., to high-T 
expansions) [l2j . Conversely, quantum fluctuations in 
the metallic Fermi-liquid regime can significantly increase 
D for T 0; this manifestation of the Kondo effect 
has been discussed in the context of Pomeranchuk cool- 
ing (without taking the low-T AF phase into account) 
[lOl l26l [27| . In contrast, the enhancement seen in Fig. 
[1] is a strong- coupling effect close to half filling as illus- 
trated in Fig.[2j In a fully developed AF state (a), a cen- 
tral spin-up atom (central black arrow) can hop virtually 
to all Z = 6 next neighbors, lowering its energy (in 2"*^ 
order) to i?AF — —Zt^/U. In contrast, half of the neigh- 
boring sites are (on average) forbidden by the Pauli prin- 
ciple in a paramagnetic state (b), thus Ep = ~Zi^ / {2U). 
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FIG. 2: (Color online) Illustration of mechanism for enhanced 
double occupancy (at strong coupling) in the AF state (a) 
compared to a paramagnetic state (b); see text, (c) DMFT- 
QMC estimates of double occupancy at half filling versus T 
for various interactions U; arrows indicate corresponding Neel 
temperatures. Thin lines: results for metastable paramag- 
netic phase, (d) same data scaled to values of critical point. 



By D ~ dE/dU (valid at T = 0), the argument implies 

Dxp/D-p 2. Thus, nearest-neighbor AF correla- 

tions can double D compared to the paramagnet at large 
U and low T . Ferromagnetic correlations, on the other 
hand, would suppress D. 

The DMFT data shown in Fig.[2](c) for a half-filled ho- 
mogeneous system confirm this (dimension independent) 
picture: At [/ = 18t (dashed line), the double occupancy 
increases strongly below T/v ~ O.BBt (arrow); the relative 
enhancement at T « 0.2t exceeds 70%, as best seen in 
the scaled view of Fig. [2] (d) . In contrast, no temperature 
dependence is visible in the paramagnetic phase, which 
may be continued artificially (thin lines) down to T = 0. 
Although the relative AF enhancement of D is smaller 
at [/ = 12i (solid lines), both the larger absolute scale 
and the larger Tn ~ 0.46 make this parameter more fa- 
vorable for experiment; whereas thermal fluctuations are 
still small at T sa Tn . The strong-coupling effect vanishes 
rapidly for smaller U\ sAU ~9t (dash-dotted line), the 
AF correlations even lead to a reduction of D. 

For quantitative predictions, let us now return to 
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FIG. 3: (Color online) Radial dependence of staggered mag- 
netization (top row) and double occupancy (bottom row) 
for V = 0.05t at temperature T = 0.25f (left column) and 
T = 0.42t (right column). The central row shows the relative 
enhancement of D in the AF phase. Thick lines correspond 
to RDMFT data; thin lines represent LDA results in (a)-(d) 
and paramagnetic calculations in (e)-(f). 



RDMFT results for inhomogeneous systems. As seen in 
Fig- El (a), the order parameter (thick lines) is nearly con- 
stant in the AF core at T = 0.25t, with mf^ll > 0.8 for 
U/t = 9, 10, 12 and m^|^ w 0.7 for U 18t, before it 
slowly decays to zero. The width of the transition region 
(3 to 4 lattice spacings) must be attributed to proximity 
effects [l6| since the phase boundaries are sharp within 
LDA (thin lines). The double occupancies (thick lines 
in panel e) show strong enhancement throughout the AF 
core in comparison with (enforced) paramagnetic solu- 
tions (thin hues) for U/t 12, 18; for U < lOt, the lat- 
ter yield even larger D in the center. As seen in Fig. [3] 
(c), the relative enhancement of D is almost equal for 
U = 12t and U = 18t at this fixed temperature (cf. Fig. 

the strong deviations between RDMFT and LDA 
solutions indicate significant proximity effects for all U. 
At larger T = 0A2t, the AF order is generally weakened 
(and disappears for U = 18t) as shown in Fig. [3] (b); the 
transitions also become smoother, especially towards the 
shrinking core. Both effects result in a smaller enhance- 
ment of D for U = 12t (Fig.lHd) and (f)). 

Let us, finally, discuss whether the double occupancy 
retains its clear signals for AF correlations in averages 
over large parts of the cloud 0] . This is indeed the case, 
as shown in Fig. [3] for U = 12t: although the total frac- 
tion i?frac = 2^jiI'i/-/V of atoms on doubly occupied 
sites (circles) is temperature dependent also at elevated 
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FIG. 4: (Color online) Fraction of atoms on doubly occu- 
pied lattice sites as a function of squared temperature; circles 
represent RDMFT results for the full cloud; polygons corre- 
spond to measurements along a gaussian beam of radius R, 
as illustrated in the inset. Thick lines are quadratic fits for 
(T/t)^ > {Tn/tf « 0.21 (and guides to the eye below); their 
extrapolations T — >■ are shown as thin lines. 

temperatures, due to the impact of the metallic shell, 
the strong increase in the AF phase at (T/t)'^ < 0.2 can 
clearly be distinguished from the extrapolated param- 
agnetic behavior (thin dotted line). As expected, the 
signal becomes sharper when the measurements are con- 
centrated [1^ towards the AF core; for R = 2a (triangles 
in Fig.ll]) the slope of the curve is negligible above Tn, so 
that the increase of -Dfrac with the onset of AF is almost 
as pronounced as in the homogeneous case (Fig. [2]). 



Conclusion - In this Letter, we have presented the first 
simulation of the Neel transition in fermions confined in 
an optical lattice at full DMFT accuracy, using a new 
QMC based RDMFT implementation for realistic trap 
parameters. We have found that the onset of AF order 
at low T is signaled in the strong-coupling regime {U > 
W = 12t) by an enhanced double occupancy. This is a 
very specific signature of antiferromagnetism since D is 
(nearly) temperature independent in the paramagnetic 
phase in this regime (for T <t/2 and within DMFT). 

It is clear that nonlocal correlations (beyond DMFT) 
will affect some of the results of this study: e.g., long- 
range AF order can be expected in 3 dimensions only 
about 25% below T^^^'^ On the other hand, DMFT 
can reproduce recent QMC measurements of D{T) in 2 
dimensions ^ (for U/t = 8^ and low T) within 10%; 
even the temperature scale Tgpin for short-range AF cor- 
relations m is only 10% below T^^^^ (at U/t = 8). 
We may therefore expect our predictions for D{T) to be 
quantitatively accurate in 3 dimensions, up to a shift of 
the kink from T^^^""" to the true Tn- In particular, the 
low-temperature enhancements of D should set in around 
j-DMFT^ as a signature of increasing short-range AF cor- 
relations. High-precision measurements of the double 
occupancy will thus be able to reveal AF correlations 
(as a precursor for AF order) at essentially the tempera- 
tures/entropies reached in current experiments. 
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